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Ice XV: a new thermodynamically stable phase of ice 
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A new phase of ice, named ice XV, has been identified and its structure determined by neutron 
diffraction. Ice XV is the hydrogen- ordered counterpart of ice VI and is thermodynamically stable 
at temperatures below ~130 K in the 0.8 to 1.5 GPa pressure range. The regions of stability 
in the medium pressure range of the phase diagram have thus been finally mapped, with only 
hydrogen-ordered phases stable at K. The ordered ice XV structure is antiferroelectric (PI), in 
clear disagreement with recent theoretical calculations predicting ferroelectric ordering (Cc). 

PACS numbers: 64.60.Cn, 61.05.fm, 64.70.kt, 61.50.Ks 
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When water freezes, hydrogen-disordered phases of ice 
crystallize which exhibit orientational disorder of the hy- 
drogen bonded water molecules [j], Upon isobaric 
cooling, these phases are expected to undergo transitions 
to hydrogen-ordered phases in which the water molecules 
adopt the energetically most favored orientations (cf. 
Fig. 1). However, because of the highly cooperative na- 
ture of molecular reorientation in ice, the sluggish kinet- 
ics of molecular reorientation at low temperatures often 
prevent these ordering transitions occurring, and the for- 
mation of 'orientational glasses' is observed instead 

The region in the phase diagram at temperatures be- 
low the stability domain of ice VI is still 'uncharted 
territory' (cf. Fig. 1(a)); it is unknown which phase 
is thermodynamically stable under these conditions 0. 
Ice VI is hydrogen-disordered, and extending its region 
of stability down to K would result in the problem- 
atic situation that a phase, with configurational entropy 
greater than zero, would be stable at K. The transi- 
tion of ice VI to a hydrogen-ordered phase upon cool- 
ing would be a way out of this dilemma. Alternatively, 
it could be that the ice II/ VI and ice VI/ VIII equilib- 
rium lines meet above K at an ice II/ VI/ VIII triple 
point below which the ordered ices II and VIII would 
be the stable phases @]. Based on unpublished neutron 
diffraction data, Kamb suggested the formation of a par- 
tially hydrogen-ordered ice VI phase [5j. Also, Johari 
and Whalley reported observing a very slow transforma- 
tion upon cooling ice VI @. However, no diffraction 
data were obtained in that study. Subsequent neutron 
diffraction measurements concluded that no significant 
structural changes were observed when pure ice VI was 
cooled to low temperatures Q. In agreement with this 
conclusion, Raman spectroscopic measurements of ice VI 
at low temperatures also did not reveal the features ex- 
pected for hydrogen ordering ||, even when the sample 



was doped with potassium hydroxide M, a recipe found 
to promote hydrogen ordering of ice Ih [1(| • Recently, it 
was shown that doping the hydrogen-disordered ices V 
and XII with minute amounts of hydrochloric acid helps 
maintaining dynamic states to low enough temperatures 
so that the phase transitions to the hydrogen-ordered ices 
XIII and XIV could be observed 
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The ordering processes in ice lead to either ferroelec- 
tric structures, with the water dipole moments adding 
up to yield a net moment, or antiferroelectric struc- 
tures in which the water dipole moments cancel each 
other. The existence of ferroelectric ice would have con- 
sequences for the geochemistry and mineral physics of 
icy comets, moons and planets [Hj]. So far, one ferro- 
electric ice phase has been observed experimentally (ice 
XI) [l(| HH. However, computational studies have pre- 
dicted that the hydrogen-ordered counterpart of ice VI 
should be ferroelectric as well [TJ, [l5| . 

In this Letter we show that hydrochloric acid doped ice 
VI transforms to a new, hydrogen-ordered phase of ice, 
named ice XV, upon cooling. The crystal structure of ice 
XV is determined by neutron diffraction and compared 
to the computational predictions [3, EH- Finally, the 
thermodynamic stability of ice XV is analysed, and the 
revised phase diagram of water and ice is presented. 

DC1 doped D2O ice VI samples were prepared by freez- 
ing a 0.01 mol L _1 solution and transferring the finely 
ground sample into a Paris-Edinburgh pressure cell |16( 
which was precooled to liquid nitrogen temperature, then 
pressurized and heated to 250 K at about 0.9 GPa. The 
crystallization of ice VI was followed in-situ by record- 
ing powder neutron diffraction patterns on the PEARL 
beamline at ISIS. The sample was then slow-cooled to 80 
K at 0.2 K min -1 under pressure and the pressure sub- 
sequently released (cf. Fig. 1(a)). After this, the sample 
was recovered from the pressure cell for further analysis 
on the GEM diffractometer. 

Figure 2(a) shows the changes of the lattice constant c 
upon heating the sample at ambient pressure from 80 K 
to 138 K. The lattice constant starts to contract at about 
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FIG. 1: (Color online) (a) The phase diagram of water and 
ice, including the liquidus lines of metastable ices IV and 
XII (long dashed lines) and extrapolated equilibria lines at 
low temperatures (short dashed lines). Hydrogen-disordered 
and ordered phases are indicated by italic and regular Roman 
numerals, respectively. DC1 doped ice VI was produced by 
heating doped ice lh to 250 K at ~0.9 GPa. The consecutive 
p/T pathways are shown, (b) Schematic illustration of the 
hydrogen ordering process as followed by diffraction. Larger 
spheres indicate oxygen atoms and smaller spheres hydrogen 
atoms, respectively. Fully occupied atomic sites are shown in 
black and half occupied sites in gray. 
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124 K upon first heating. This change was found to be re- 
versible upon slow-cooling to 80 K at 0.1 K min -1 , with c 
expanding discontinuously starting at about 130 K. Sub- 
sequent heating showed again reversible behavior, and c 
reached the value obtained after the first heating. The 
stepwise change upon second heating was almost three 
times larger than observed upon first heating. It was re- 
cently shown that hydrogen order / disorder phase tran- 
sitions in ices may be accompanied by such discontinu- 
ous changes in the lattice constants and that the magni- 
tude of a step correlates with the degree of order in the 
hydrogen-ordered phase [ll|, [TtJ ■ The behaviour shown 
in Fig. 2(a) therefore suggests that a hydrogen order / 
disorder phase transition may be occurring. 

The inset in Fig. 2(b) shows the powder diffraction 
data recorded between the various heating and cooling 
steps. For the sample after slow-cooling under pressure, 
Bragg peaks (indicated by arrows) were observed which 
are either not permitted by the space group symmetry 
of ice VI, P42 1 rime, or which have zero intensity in that 
structure. These peaks disappeared upon heating to 138 



FIG. 2: (Color online) Analysis of the powder neutron diffrac- 
tion data (GEM instrument), (a) Changes of lattice constant 
c upon heating and cooling at ambient pressure. The num- 
bers in (a) and (b) correspond to the p/T waypoints in Fig. 
1(a). (b) Observed, calculated and difference neutron powder 
diffraction profiles for ice XV after slow-cooling to 80 K at 
ambient pressure (R wp = 0.0155, R p = 0.0186). The crystal 
structure was refined using the GSAS program fl8j ] and space 
group PI. Linear constraints and chemical composition re- 
straints were used to maintain the stoichiometries of the hy- 
drogen bonds and water molecules, respectively. Tick marks 
indicate the positions of the Bragg peaks for ice XV, Ni and 
V (top to bottom), (inset) Magnification of a region where 
new Bragg peaks appeared. The data shown were recorded 
at 80 K after slow-cooling under pressure (bottom), at 138 K 
after heating at ambient pressure (middle) and at 80 K af- 
ter slow-cooling at ambient pressure (top). Arrows indicate 
the positions of new Bragg peaks, (c) \ 2 values obtained after 
converged Rietveld refinements using the fully ordered models 
shown in Fig. 3(b) and space group symmetries Cc, Pn, P2i 
and Pi. The \ 2 value obtained for the structure predicted in 
refs [111 , fist ] is depicted as an open square. 



K and reappeared upon slow-cooling to 80 K. In fact, af- 
ter slow-cooling at ambient pressure, the new peaks were 
more intense than after slow-cooling under pressure. In 
agreement with the changes observed of the lattice con- 
stant c (c/. Fig. 2(a)), this further suggests the forma- 
tion of a more highly hydrogen-ordered structure after 
slow-cooling at ambient pressure. 

The tetragonal unit cell of ice VI comprises ten water 
molecules belonging to two independent but interpene- 
trating hydrogen bonded networks (c/. Fig. 3(a)) Q. 
Hence, the ice VI structure can be regarded as a 'self- 
clathrate' @, ■ The central structural motifs of each 
network are 'water hexamers' which share corners along 
the crystallographic c direction and are hydrogen bonded 
to each other along the a and b directions. The P4 2 /nmc 
space group symmetry requires ice VI to be fully hydro- 
gen disordered, which means that the fractional occupan- 
cies of the hydrogen positions are restricted to \ (7L Il9l] . 
Consequently, any change in hydrogen order has to be 
associated with breaking of the space group symmetry 
i. 

For a single hydrogen bonded network there are six 
different ways in which the structure can become hydro- 
gen ordered as shown in Fig. 3(b). Due to the chiral- 
ity of the ordered hexameric units, the structures can be 
grouped into pairs which are related by mirror symmetry 
(A1/A2, B1/B2 and C1/C2). The A1/A2 and B1/B2 
structures are polarized in the a-b plane only, whereas 
the C1/C2 structures show additional polarization in the 
c direction. The space groups which can accommodate 
two symmetry-related, fully ordered networks are Cc, Pn, 
P2i and PI. Complete cancelation of the polarities of the 
two networks, and thus antiferroelectrity, is only achieved 
in PI. The other space group symmetries yield ferroelec- 
tric structures. Space group Cc can be excluded conclu- 
sively due to the presence of a new Bragg peak at 1.93 A 
which is forbidden in Cc (cf. inset in Fig. 2(b)). Simi- 
larly, the new weak peak at 2.08 A cannot be indexed in 
the Pn space group. 

Our powder diffraction data collected after slow- 
cooling at ambient pressure were initially fitted by per- 
forming Rietveld refinements using the different possible 
fully ordered structural models. The \ 2 values obtained 
after convergence of the refinements indicate the quality 
of the fits and are shown in Fig. 2(c). Consistent with the 
absence of the new Bragg peaks, the structural models 
with Cc and Pn space group symmetries produce unsat- 
isfactory fits. For the P2i and Pi models, a clear prefer- 
ence emerges for ordered structures containing the CI or 
C2 enantiomers, either as a true racemate (PI) or as a 
racemic mixture (P2i). For the fully ordered structures, 
the best fits in P2i and PI are of similar quality; how- 
ever, PI yields a significantly better agreement with the 
data when the occupancies of the hydrogen atom sites are 
refined. For the P2i structures, the occupancies of some 
of the hydrogen positions increased above 1 or fell below 
during the refinement while others converged close to 
i, behavior which indicates incorrect structural models. 
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(c) ice XV 
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FIG. 3: (Color online) Unit cell projections of (a) hydrogen- 
disordered ice VI (P42/nmc), (b) the six possible hydrogen- 
ordered structures of a single hydrogen bonded network (PI) 
and (c) hydrogen-ordered ice XV (PI). The atomic site occu- 
pancies are indicated by a grayscale gradient from black (fully 
occupied) to white (not occupied), and the polyhedra show 
the front faces of the hexameric units. The oxygen atoms be- 
longing to different networks are marked with '1' and '2' in 
(a). 



We therefore assign PI space group symmetry to the new 
hydrogen-ordered phase. Unlike P2i , space group PI can 
be obtained by a single irreducible representation of the 
high-symmetry space group P^jnmc. Therefore, Lan- 
dau theory [20|, |2l| permits a one-step, weakly first order 
transition such as the one observed here, where the or- 
der parameter is a scalar field representing the departure 
from half occupancy of the hydrogen sites. Out of all the 
possible ordered structures, the experimental structure 
contains the most polarized networks {CI and C2) and 
its space group symmetry allows for the most effective 
canceling of the network polarities. 

On the basis of the reduced space group symmetry, and 
the reversible and discontinuous change in the lattice con- 
stant, we identify this as a new phase of ice to which we 
assign the Roman numeral XV. Table I shows the refined 
atomic site occupancies, which are, consistent with the 
changes in the lattice constant, more ordered after cool- 
ing at ambient pressure than after cooling under pres- 
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TABLE I: Fractional coordinates, isotropic atomic displace- 
ment parameters and fractional occupancies for ice XV (PI). 
Data were collected at 80 K after cooling the sample at 0.1 
K min -1 from 138 to 80 K at ambient pressure. Deuterium 
positions with refined fractional occupancies < | are omitted 
from table, a = 6.2323 A, b = 6.2438 A, c = 5.7903 A, a 
= 90.06°, /3 = 89.99°, 7 = 89.92°. The occupancy values in 
parentheses are for the sample after cooling at 0.2 K min -1 
at ~0.9 GPa. 
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sure, and therefore closer to the thermodynamic limit of 
complete hydrogen order. The psewrfo-orthorhombic unit 
cell of ice XV with slightly distorted hexameric units is 
shown in Fig. 3(c). 

Density functional theory approaches have successfully 
predicted the structure of hydrogen-ordered ice XIV [22( 
and have produced results consistent with the experimen- 
tal ice XIII structure being the energetic ground state 
in the appropriate p/T ranges [23| . To our knowledge, 
there are two independent computational predictions of 
the hydrogen-ordered structure of ice XV [3, EH ■ Both 
studies propose that the hydrogen-ordered ground state 
is ferroelectric with Al / A 2 hydrogen-ordered networks 
and Cc space group symmetry, which is in clear dis- 
agreement with the experimental structure reported here. 
The reasons for this disagreement are not understood at 



present. 

Upon isobaric cooling, the increase in free energy of 
a phase correlates with its entropy, (§p) P = —S. It fol- 
lows that at temperatures below the ordering tempera- 
ture (~130 K from Fig. 2(a)), hydrogen-ordered ice XV 
will have a lower free energy than ice VI [24(. Since the 
ordering transition takes place in what was previously 
thought to be the stability region of ice VI, the lower tem- 
perature part of the ice VI region of stability is now as- 
signed to the stability region of ice XV. The phase bound- 
aries of ice XV with ices II and VIII can be constructed 
by considering the Clapeyron equation, (^)ag = 
which requires approximately vertical phase boundaries 
between two hydrogen-ordered phases of different den- 
sity (II/XV and XV/ VIII), and approximately horizontal 
phase boundaries between a hydrogen-disordered and a 
hydrogen-ordered phase of similar density (VI /XV) . The 
shaded area in Fig. 1(a) indicates the region in the phase 
diagram in which ice XV is expected to be thermodynam- 
ically more stable than any other of the currently known 
phases. 

In summary, we have prepared a new, hydrogen- 
ordered phase of ice, named ice XV, and solved its crystal 
structure. Ice XV is antiferroelectric, thermodynamically 
stable at low temperatures in the 0.8 to 1.5 GPa pres- 
sure range and forms upon slow-cooling hydrochloric acid 
doped ice VI. The revised phase diagram is now finally 
consistent with only hydrogen-ordered phases of ice being 
the stable phases at K. The crystal structure of ice XV 
once again underlines the importance of experimental re- 
sults to serve as benchmarks for computational models of 
water and to test their abilities to reproduce the lowest 
energy states. Such studies could lead to improved theo- 
retical descriptions of the interactions of water with itself, 
and with more complex inorganic and organic molecules. 
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(C.G.S.). 
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